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Research  Accomplishments 

This  project  was  funded  in  March  1,  1995  for  work  on  three  problems  in  two-phase  flow. 

Problem  One 

The  first  problem  was  to  simulate  the  motions  of  particles  in  a  viscoelastic  fluid  in  a  direct,  two- 

dimensional  simulation.  The  main  results  are  described  in  the  following  three  abstracts. 

Dynamic  Simulation  of  sedimentation  of  solid  particies  in  an  Oidroyd-B  fiuid 

/.  Fen^  P.Y.  Huang,  D.D.  Jos^h 

Abstract 

In  this  paper  we  present  a  two-dimensional  numerical  study  of  the  viscoelastic 
effects  on  the  sedimentation  of  particles  in  the  presence  of  solid  walls  or  another 
particle.  The  Navier-Stokes  equations  coupled  with  an  Oldroyd-B  model  are 
solved  using  a  finite-element  method  with  the  EVSS  formalism,  and  the  particles 
are  moved  according  to  their  equations  of  motion.  In  a  vertical  channel  filled 
with  a  viscoelastic  fluid,  a  particle  settling  very  close  to  one  side  wall  experiences 
a  repulsion  from  the  wall;  a  particle  farther  away  from  the  wall  is  attracted 
toward  it.  Thus  a  settling  particle  will  approach  an  eccentric  equilibrium  position, 
which  depends  on  the  Reynolds  an  Deborah  numbers.  Two  particles  settling  one 
on  top  of  the  other  attract  and  from  a  doublet  if  their  initial  separation  is  not  too 
large.  Two  particles  settling  side  by  side  approach  each  other  and  the  doublet 
also  rotates  till  the  line  of  centers  is  aligned  with  the  direction  of  sedimentation. 

The  particle-particle  interactions  are  in  qualitative  agreement  with  experimental 
observations,  while  the  wall  repulsion  has  not  been  documented  in  experiments. 

The  driving  force  for  lateral  migrations  is  shown  to  correlate  with  the  pressure 
distribution  on  the  particle’s  surface.  As  a  rule,  viscoelasticity  affects  the  motion 
of  particles  by  modifying  the  pressure  distribution  on  their  surface.  ITie  direct 
contribution  of  viscoelastic  normal  stresses  to  the  force  and  torque  is  not 
important. 
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Direct  simulation  of  the  motion  of  solid  particles  in  Couette  and  Poiseuille  flows 
of  viscoelastic  fluids 

P,Y,  Huang,  J.  Feng  HH,  Hu  and D.D,  Joseph 
Abstract 

This  paper  reports  the  results  of  direct  numerical  simulation  of  the  motion  of  a 
two-dimensional  circular  cylinder  in  Couette  flow  and  Poiseuille  flow  of  an 
Oldroyd-B  fluid.  Both  neutrally  buoyant  and  non-neutrally  buoyant  cylinders  are 
considered.  The  cylinder's  motion  and  the  mechanisms  which  cause  the  cylinders 
to  migrate  are  studied.  The  stable  equilibrium  position  of  neutrally  buoyant 
particles  varies  with  inertia,  elasticity,  shear-thinning  and  the  blockage  ratio  of 
the  channel  in  both  shear  flows.  Shear-thinning  promotes  the  migration  of  the 
cylinder  to  the  wall  while  inertia  causes  the  cylinder  to  migrate  away  from  the 
wall.  The  cylinder  moves  closer  to  the  wall  in  a  narrower  channel.  In  a  Poiseuille 
flow,  the  effect  of  elastic  normal  stresses  is  manifested  by  an  attraction  toward 
the  nearby  wall  if  the  blockage  is  strong.  If  the  blockage  is  weak,  the  normal 
stresses  act  through  the  curvature  of  the  inflow  velocity  profile  and  generate  a 
lateral  force  that  points  to  the  centerline.  In  both  cases,  the  migration  of 
particles  is  controlled  by  elastic  normal  stresses  which  in  the  limit  of  slow  flow 
in  two  dimensions  ate  compressive  and  proportional  to  the  square  of  the  shear 
rate  on  the  body.  A  slightly  buoyant  cylinder  in  Couette  flow  migrates  to  an 
equilibrixam  position  nearer  the  centerline  of  the  channel  in  a  viscoelastic  fluid 
than  in  a  Newtonian  fluid;  On  the  other  hand,  the  same  slightly  buoyant  cylinder 
in  Poiseuille  flow  moves  to  a  stable  position  farther  away  from  the  centerline  of 
the  channel  in  a  viscoelastic  fluid  than  in  a  Newtonian  fluid.  Marked  effects  of 
shear  thinning  are  documented  and  discussed. 


Direct  simulation  of  the  sedimentation  of  elliptic  particles  in  Oldroyd-B  fluids 

■P.  Y.  Huang,  H.H.  Hu  and  D.D.  Joseph 

Abstract 

Cross  stream  migration  and  stable  orientations  of  elliptic  particles  falling  in  an 
Oldroyd-B  fluid  in  a  channel  are  studied.  We  show  that  the  normal  component 
of  the  extra  stress  on  a  rigid  body  vanishes;  lateral  forces  and  torques  are 
determined  by  the  pressure.  Inertia  turns  the  longside  of  the  ellipse  across  the 
stream  and  elasticity  turns  it  along  the  stream;  tilted  off-center  falling  is  unstable. 

There  are  two  critical  numbers;  elasticity  and  Mach  numbers.  When  the  elasticity 
number  is  smaller  than  critical  the  fluid  is  essentially  Newtonian  with  broadside- 
on  falling  at  the  centerline  of  the  channel.  For  larger  elasticity  numbers  the 
settling  turns  the  longside  of  the  particle  along  the  stream  in  the  channel  center 
for  all  velocities  below  a  critical  one,  identified  with  a  critical  Mach  number  of 
order  one.  For  larger  Mach  numbers  the  ellipse  flips  into  broadside-on  falling 
again.  The  critical  numbers  are  functions  of  the  channel  blockage  ratio,  the 
particle  aspect  ratio  and  the  retardation/relaxation  time  ratio  of  the  fluid.  Two 
ellipses  falling  nearby,  attract,  line-up  and  straighten-out  in  a  long  chain  of 
ellipses  with  longside  vertical,  all  in  a  row.  Stable,  off-center  tilting  is  found  for 
ellipses  falling  in  shear  thinning  fluids  and  for  cylinders  with  flat  ends  in  which 
particles  tend  align  their  longest  diameter  with  gravity. 

Problem  Two 

The  second  problem  was  to  identify  new  physics  that  are  revealed  about  the  motion  and  mixing 
when  account  is  taken  of  the  fact  that  the  density  of  mixtures  of  incompressible  liquids  must  change 
when  they  are  diffusing.  These  ideas  are  thoroughly  explored  in  the  abstract  of  the  following  paper. 

Non-solenoidal  velocity  effects  and  Korteweg  stresses  in  simple  mixtures  of 
incompressible  liquids 

Daniel  D.  Joseph,  Adam  Huangs  Howard  Hu 

Abstract 

We  study  some  basic  problems  of  fluid  dynamics  of  two  incompressible  miscible 
liquids  modeled  as  a  simple  mixture  in  which  the  volume  of  the  mixture  does 
not  change  on  mixing.  In  general,  the  expansion  A  =  divu  in  these  problems 
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does  not  vanish.  The  velocity  in  such  a  mixture  can  be  decomposed  into  a 
solenoidal  and  an  expansion  part.  The  expansion  velocity  is  induced  by  diffusion 
which  is  proportional  to  the  gradient  of  the  volume  fraction  in  a  simple  mixture. 

The  expansion  can  be  large  at  certain  times  and  places.  We  have  carried  out  an 
analysis  of  transient  or  dynamic  interfadal  tension  for  the  problem  of  smoothing 
of  an  initial  discontinxiity  of  composition  across  a  plane  and  spherical  surface. 

The  dynamic  tension  at  the  spherical  interface  decays  as  it  has  two  terms, 
one  term  arises  from  the  Korteweg  stress  and  it  gives  rise  to  a  stress  opposing 
the  internal  pressure  as  in  the  case  of  equilibrium  tension  if  the  Korteweg 
coefficient  has  the  appropriate  sign.  The  other  term  arises  from  the  expansion 
velocity  and  is  proportional  to  the  rate  of  change  of  viscosity  with  volume 
fraction.  This  term  has  the  wrong  sign  for  interfacial  tension  in  the  case  of 
glycerin  and  water  solutions  but  has  the  right  sign  when  the  light  fluid  is  more 
viscous.  In  the  context  of  the  new  theory,  we  derive  a  new  and  elementary 
solution  which  describes  diffusion  of  binary  species  along  a  pipe.  An  analysis  of 
this  solution  to  small  disturbances  leads  to  a  nonseparable  problem  which  is 
solved  by  a  finite  element  method.  The  numerical  study  indicates  stability  under 
all  circumstances. 

We  found  that  diffusion  layers  ‘develop  at  side  walls  when  the  densities  of  incompressible 

constituents  of  miscible  liquids  are  different.  This  is  explained  in  the  abstract  of  the  following  paper. 

Sidewall  effects  in  the  smoothing  of  an  initial  discontinuity  of  concentration 

T.Y.  Uao  and D.D.  Joseph 

Abstract 

The  velocity  field  of  a  binary  mixture  of  incompressible  miscible  liquids  is  non- 
solenoidal  when  the  densities  of  the  two  liquids  are  different.  If  the  mixture 
density  is  linear  in  the  volume  fraction,  as  in  the  case  of  simple  (ideal)  mixtures 
or  very  nearly  for  glycerin  and  water,  then  the  velocity  can  be  decomposed  into  a 
solenoidal  and  an  expansion  part.  In  the  context  of  this  theory,  we  derive  a  new 
solution  which  describes  the  smoothing  of  an  initial  plane  discontinuity  in 
concentration  across  a  channel  bounded  by  side  walls.  The  requirement  that  the 
velocity  vanishes  on  the  side  wall  introduces  a  different  initial  discontinuity  not 
present  in  the  solenoidal  theory.  The  problem  may  be  reduced  to  a  partial 


7 


differential  equation  in  two  similarity  variables  one  for  the  smoothing  of  a 
concentration  discontinuity  without  sidewalls  and  the  other  for  the  smoothing 
the  velocity  discontinuity  at  the  sidewall.  The  similarity  equations  are  solved 
explicitly  in  a  special  case. 

Problem  Three 

The  third  problem  was  to  study  breakup  of  thickened  liquids  at  high  Mach  numbers.  We  built  a 
shock  tube  to  study  this  problem  and  found  funds  from  an  NSF  equipment  grant  to  buy  a  drum 
camera.  We  made  the  first  movies  of  drop  breakup  in  these  extreme  conditions.  It  appears  that  the 
shock  tube  is  the  only  way  presently  known  to  obtain  drop  breakup  data  at  low  pressures.  The 
results  are  reported  in  the  following  abstract. 

Breakup  of  a  liquid  drop  suddenly  exposed  to  a  high-speed  airstream 

Bj  Daniel  D.  ]osephy  /.  Belanger  G.S,  Beavers 

Submitted  to  /.  Fluid  Mech 

Abstract 

The  problem  is  to  study  how  drops  of  various  liquids,  Newtonian  and  non- 
Newtonian,  break  in  the  high  speed  akstream  created  by  a  shock  wave  in  a  shock 
tube.  We  constoicted  high  speed  movies  of  breakup  with  a  rotating  drum 
camera  giving  one  photograph  every  5|Xs.  Drops  of  the  order  of  one  millimeter 
are  reduced  to  droplet  clouds,  and  possibly  to  vapor,  in  a  time  less  than  SOOjas. 

The  entire  fragmentation  history  is  now  available  in  movie  form 
(http:/ /www.aem.umn.edu/research/Aerodynamic_Breakup).  The  movies 
reveal  sequences  of  breakup  events  which  have  universal  features  which  occur  in 
all  the  liquids  and  specific  features  which  differ  from  one  type  of  liquid  to 
another.  In  particular,  the  viscoelastic  fluids  give  rise  to  stringKke  and  clumped 
drop  fragments  not  seen  in  Newtonian  flmds.  Rayleigh-Taylor  instabilities  are 
endemic  to  the  breakup;  for  example,  the  acceleration  from  gas  to  liquid 
required  for  this  instability  is  of  the  order  10^  times  the  acceleration  of  gravity  in 
the  flow  behind  Mach  3  shock  waves.  The  process  of  drop  breakup  can  be 
characterized  in  part  as  the  inflation  of  the  drop  by  pumping  in  fingers  of  hot  ak 
under  the  action  of  Rayleigh-Taylor  instability. 
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